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In the present work the so-called Higher Tamm-Dancoff Apporximation method is pre- 
sented for the generalized case of isovector and isoscalar residual interactions treated 
simultaneously. The role of different particle- hole excitations and of proton-neutron pair- 
ing correlations in the ground state of the self-conjugate 64 Ge nucleus is discussed. 
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1. Introduction 

It was recognized already in the sixties that the proton-neutron pairing may play a 
non-negligible role in the nuclei with equal numbers of protons and neutrons due to 
the particulary large spatial overlaps of the wave functions of both types of nucleons. 
One of the first methods developed to take into account such correlations were BCS- 
type approaches. 1,2,3,4,5,6 It is however known that the BCS wave function violates 
the particle number conservation which is especially harmful in the cases of weak 
pairing and phase transitions. A generic feature of the results of BCS calculations 
generalized so as to accomodate neutron-proton correlations (i.e. that take into 
account isoscalar and isovector pairs) in self-conjugate nuclei is a sharp transition 
from like-particle to T = pairing modes. Such a scenario is absent in the particle 
number conserving approaches which allow for the coexistence of different superfluid 
phases. 7,8 ' 9 

In the present study we propose an approach which is free of the deficiencies 
of the treatments a la Bogoliubov. This method, called Higher Tamm-Dancoff Ap- 
proximation (HTDA), was already applied to study ground and isomeric states in 
heavier nuclei and to describe pairing correlations in high spin states. 10,11 Here we 
extend the HTDA formalism to take into account the proton-neutron pairing. 
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2. The HTDA framework 

The purpose of the method is to describe the pair diffusion phenomenon around the 
Fermi surface as a particular, yet very important, part of some many quasiparticle 
excitations of the particle- hole type over a Slater determinant vacuum l^o)- Thus 
one can consider dealing with a kind of higher Tamm-Dancoff approximation or 
equivalently with a highly truncated shell-model calculation. The latter point of 
view makes it clear that the success (i.e., the fast convergence in terms of the 
number of considered quasiparticles) of such a truncation scheme will depend on 
the realistic character of the vacuum in use. Such a favorable situation may be 
expected when the vacuum is defined self-consistently from the mean-field carrying 
most of the single-particle (sp) properties (even those yielded by the correlations) 
associated with a given effective Hamiltonian 

H = K + V, (1) 

where K is the kinetic energy and V an effective interaction. It is advisable to choose 
for the wave function the Hartree-Fock solution coming from the HF potential 
Vhf self-consistently obtained from the many-body reduced density matrix p of 
the correlated solution for the desired number of particles and possibly taking into 
account various constraints (e.g., deformation) and symmetries (e.g., time- reversal 
symmetry). We have thus 

#hf|*o) =£b|*o), (2) 

where 

H HF = K + Vhf , (3) 

Vhf being the one-body reduction of V for |^o)- 

The l^o) vacuum may now serve to construct an orthonormal N-body basis 
in which we will diagonalize Hamiltonian (1). In principle to build this basis we 
should consider the Slater determinant corresponding to the OpOh state and all the 
particle-hole excitations: lplh l^i), 2p2h ^2), 3p3h |* 3 ) and so on. It is however 
obvious that in practice we need to truncate the space of particle-hole excitations. 
Based on the former studies in the HTDA framework 10 we may assume that the 
components of the pair excitation type will be dominant in the ground state. In 
the section 3 we will show for the case of pp and nn correlations that the other 
components are indeed of a minor importance. 

The total ground state wave function can be decomposed in the following way 
(r and r' denoting two different charge states) 

I*) = |* T ®* T ') =xoo\%®%) + Xio|*i <S>*o')+ Xoi|*o®*i'> 

(lplh) T (lplh) T / 
(lplh) T (lplh) T , (2p2h) T (2p2h) T , 

+ ••• (4) 
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The ensemble of products of Slater determinants ® 'J J ) represents a com- 
plete orthogonal basis of the N neutrons and Z protons physical space with real 
coefficients Xm fulfilling the relation 

£Xm = l (5) 

m 

which guarantees the normalization of the function (4). It is clear that this function 
has a good particle number (^\ N\*S?) = A. To obtain the correlated ground state one 
should perform the diagonalization of the Hamiltonian (1) in the given many-body 
basis. 

Let us rewrite the Hamiltonian (1) in the following form 

H = K + Vuf- (*o|Vl*o> + V-Vhf + (*o|Vl*o) 

= (^ \H\^ )+H lQP + V Ies , (6) 

where the independent quasiparticle Hamiltonian Hiqp reads 

HiQ P = J2&vlvi, (7) 

i 

where rjj is the particle creation operator a\ for i being a particle (unoccupied) 
state and to the anihhilation operator en in the case of hole (occupied) states. For 
we have: £j = e l or & = — e l , with e 1 indicating the single-particle (HF) energy 
corresponding to the particle or hole case, respectively. As seen from Eq. (6), the 
residual interaction reads 

= v - vW + <*o|vr*o> ■ (8) 

The matrix element of the above given Hamiltonian in the multi-particle multi- 
hole basis takes the form 

Hij = ( <*o|#|*0> + £ E?- b ) S tJ + (^l^esl^) , (9) 
\ r—p^n / 

where ^2 T=p n E^_ h is the total particle-hole excitation energy calculated with re- 
spect to the vacuum |* ) = l^o ® *o )■ The residual interaction comprising only 
one- and two-body operators will merely contribute to the matrix elements in cases 
where the two many-body states are identical or differ by one or two particles. 
The many-body matrix elements of V rcs can be evaluated with the use of the Wick 
theorem (for details see Ref. 12 ) and expressed as sums of proper two-body matrix 
elements of the interaction V. 

In actual calculations with our approach we have used the Skyrme force for 
the HF part of the problem. Since the majority of the Skyrme forces is not well 
suited to reproduce the data in the particle-particle channel we apply a (5-force to 
define the residual interaction. Using the zero-range force is well established in the 
context of like-particle and- to a lesser extent- of the proton-neutron pairing. To 
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assess the validity of our approach at this early stage of the investigation of the 
proton-neutron correlations within the HTDA framework a substitution 

V=>V S = V (T)S(r 12 )n s n T (10) 

T=0,1 

is therefore done for the residual interaction. The operators n s n T project onto 
spin-isospin subspaces 

fl s = I(l-(-l)Sp-), (ll) 

n T = l(i -(-if PI, (12) 

where P a = 1/2(1 — <?i • (T 2 ), P t = 1/2(1 — fi - 7^) are the standard spin and isospin 
exchange operators. 

3. Results 

The main challenge in the exact pairing diagonalizations methods and shell-model 
like calculations is the large dimension of the matrix to be diagonalized. In the 
present study we use the Lanczos algorithm 13 within the code developed by Parlett 
and Scott 14 which allows to search for the lowest energy solutions in a reasonable 
time. However, as will be shown in the following, the advantage of the HTDA 
method is that the amount of the many-body configurations within a given single- 
particle space that need to be taken into account is quite limited and the matrix 
diagonalization problem is far from attaining the difficulty level of typical shell- 
model calculations. 

The calculations within the HTDA framework which we will discuss here are per- 
formed for the ground state of the 64 Ge nucleus as obtained in a Skyrme-HF+BCS 
approach. We use the SIII force for the particle-hole channel, seniority force in 
the BCS part and volume 6 interaction to calculate two-body matrix elements in 
the HTDA part. The HTDA results presented here are non self-consistent. Indeed, 
only one diagonalization of the HTDA Hamiltonian matrix has been carried out 
for a Hamiltonian Hhf obtained with the density p corresponding to the HF+BCS 
calculation. 

3.1. Limiting case: pp and nn interactions 

In the limiting case of vanishing proton-neutron correlations we obtain two decou- 
pled problems, one for each charge state. Let us discuss first the aspects of the basis 
truncation, both in single-particle and many-particle many-hole excitations terms. 
If the assumptions of the BCS theory are to hold, we may expect that the most 
important part of the ground state correlations will correspond to the pair exci- 
tations, i.e. the situation where two particles in Kramers degenerate hole orbitals 
are promoted to two Kramers degenerate particle levels. In what follows we will 
restrict the many-body states first to take into account all possible lplh and 2p2h 
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Table 1. Numbers of many-body configurations of different types depending on the number of sp 
levels in the window. 



Number of sp levels lplh 1 pair 2p2h 2pairs 3 pairs 



12 


10 


40 


260 


280 


560 


20 


22 


96 


1744 


1848 


12157 


28 


48 


180 


6596 


6758 


58404 


36 


76 


260 


13706 


14230 


>100000 



excitations including thus the one pair excitations in addition to the particle-hole 
vacuum. In further calculations only the states corresponding to one, two and three 
pairs excitations will be considered. 

Having chosen a configuration space size in terms of the complexity of many- 
particle many-hole states, one has to further truncate on single-particle levels from 
which these configurations will be built. We are facing thus a situation met in 
customary BCS calculations. Typically, we will limit our sp subspace to the so- 
called configuration space window defined as 

e F ±EZ, (13) 

where ep is the Fermi energy (defined as the average between the single-particle 
energies of the last occupied and first unoccupied levels) and £^ t the cut-off energy 
The actual value of the latter should be chosen such that the inclusion of additional 
sp levels does not introduce any significant physical consequences but only a possible 
slight renormalization of the relevant quantities. In Table 1, to exemplify the matrix 
dimensions, we give numbers of configurations corresponding to different particle- 
hole excitations depending on the number of sp levels included in the configuration 
space window. 

Since in the HTDA calculations we are using a 5-forcc it is clear that both 
the coupling constant and the cut-off parameter are necessary to define fully the 
interaction. It is then indispensable to fix the 5-force strength in our approach 
to reproduce physical quantities, e.g., the phenomenological 3-point gaps for the 
considered nucleus (A„ = 1.48 MeV, A p =1.14 MeV evaluated with experimental 
masses of Rcf. 15 ). Assuming that the appearance of the pairing gap is related to a 
breaking of the Cooper pair of the lowest energy, we perform the HTDA calculations 
with one level (neutron or proton) closest to the Fermi energy blocked and adopt 
the difference of the expectation value of V les in the two calculations as a proper 
measure of the pairing correlations that can be compared to the experimental odd- 
even mass differences. Namely, we define 



A = [E{n) - E IQP (n)} - [E(n - 1) - E 1QP (n - 1)] , (14) 
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Table 2. Percentage of different components of the ground state correlated wave function. The 
results concern the case where all types of particle-hole excitations up to the 2 nd order are taken 
into account. 





OpOh 


lplh 


1 pair 


2p2h (all included) 


neutrons 


69.6 


0.003 


29.6 


30.4 


protons 


70.6 


0.005 


29.3 


28.3 



where E = {9\H\^>), -Eiqp = (\f , |f/iQp|\f') and n indicates the number of sp levels 
in the configuartion space window. Since we deal with the N = Z nucleus one value 
of the coupling constant Vo and of the cut-off energy E^ t = lOMeV is adopted for 
both charge states in further calculations. 

First, let us discuss the result of the diagonalization of the Hamiltonian matrix in 
the many body basis containing all lplh and 2p2h components. The decomposition 
of the correlated function for neutrons and protons is given in Table 2. As expected, 
the main contribution beyond the OpOh vacuum comes from one pair excitations. 
The lplh part is fully negligible while the 2p2h excitations of other types constitute 
less than 1% of the correlated wave function so restricting the many-particle many- 
hole excitations to study only pair excitations in addition to the vacuum seems 
well justified. Of course it allows for a great simplification from the calculational 
point of view, as the number of pair excitations is very modest as compared to the 
number of all possible particle-hole excitations (cf. Tab. 1). 

The same procedure is employed in the case of calculations where one, two and 
three pairs excitations in addition to the vacuum are considered. In Table 3 the 
decomposition of the correlated wave function is given for neutrons and protons in 
three cases: when only 1 pair excitations are considered in addition to the vacuum 
state, when 2 pairs are added and finally, when 3 pairs are as well embedded in 
the calculations. We also give the corresponding values of the correlation energy 
defined as the difference of the mean values of the Hamilton operator of the system 
in correlated and uncorrelated states 

£corr = (*|#|*> - (*o|#|*o) ■ (15) 

As seen, one pair excitations contribute the most in all cases, two pairs excita- 
tions give only about 2% of the total amount while the probability for three pairs 
is negligible. Including the latter has a minimal influence on both the correlation 
energy and on the percentage of other components. However, the presence of two 
pairs boosts the population of the one pair excitation content of the correlated 
wave function, as already observed in the calculations of Ref. 12 . This entails that 
to describe accurately the one pair excitation phenomenon two pair transfers are a 
priori needed. Nevertheless, at this early stage of our study of proton-neutron pair- 
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Table 3. Percentage of various components of the correlated wave function and the correlation 
energy values for neutrons and protons. Results of three calculations are reported: with 1 pair 
excitations added to the vacuum, with 1 and 2 pairs and with 1, 2 and 3 pairs included. 



OpOh 1 pair 2 pairs 3 pairs E COII (MeV) 



neutrons 65.4 34.6 -1.16 
49.55 48.5 1.95 -1.43 
49.0 48.9 2.1 0.02 -1.45 



protons 65.3 34.7 -1.18 
47.7 50.2 2.1 -1.48 
46.9 50.7 2.3 0.04 -1.50 



ing correlations we will exclude 2 pairs components for obvious practical reasons. 



3.2. General case: T = 1 and T — pairing interaction 

As we already stated we aim at including the T z = pairing correlations within 
the HTDA framework. As above discussed in this preliminary study we limit the 
configuration space to the so-called one pair excitations in all channels. 

Since proton-neutron pairs may exist in both T = and T = 1 channels the 
problem of adjusting the strength of the isoscalar pairing arises here. The interplay 
of T = and T = 1 pairing being rather unclear, we will adopt the same value of 
the isovector pairing strength for pp, nn and pn channels and will treat the T = 
interaction strength as a free parameter. a 

Within a given configuration space and T — 1 strength adjusted as described 
in the precedent we perform the calculations as a function of the increasing T = 
coupling constant. In Fig. 1 the evolution of the probability of different components 
is plotted against the ratio 

x = V (T = 0)/V (T=l). (16) 

The probability of the T z — pair excitation becomes siginificant with respect to 
that of pp or nn pair excitation for x > 2 which corresponds to the ratio value 
at which a pn superfluid solution is found in BCS(<5)-type of calculations in this 
nucleus. 9 The gain in energy due to the pn coupling is non-negligible which can be 

a A considerable amount of pn pairing correlations can be as well obtained in an isospin broken 
model with VJ" T > V TT , nevertheless we do not see any argument supporting such a choice. 
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vacuum 
pn 
nn 

PP 



12 3 4 

x 

Fig. 1. Total probability of different type of components (vacuum- (0p0h) p ®(0p0h)n, pn- 
(lplh) p ®(lplh)„, pp-(0p0h) n ®(2p2h) p , nn-(2p2h) n (g)(0p0h) n ) in the correlated wave function 
versus the ratio x of the 5-force strengths in isoscalar (T = 0) and isovector (T = 1) channels. 

seen from Fig. 2 where the contribution of different components to the correlation 
energy is depicted. It is worth noting that the increase of the energy is of a similar 
magnitude like that of BCS calculations with the approximate particle number 
projection of the Lipkin-Nogami (LN) 9 nonetheless a certain gain in the HTDA 
correlation energy due to the pn coupling occurs already in the range of parameters 
for which no pn pairing solution is found in the BCS+LN method (x < 2). 

4. Summary 

We have presented a method explicitely conserving particle number in the context 
of describing isovector and isoscalar pairing correlations. The inclusion of T = 
pairs leads to a constant increase of the absolute value of the correlation energy. 
The expected benefit of the HTDA framework in yielding a much smoother decrease 
of correlations than when using a Bogoliubov quasiparticlc vacuum in weak pairing 
regimes is confirmed. 
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Fig. 2. Absolute value of the correlation energy plotted against the ratio x of the interaction 
strengths in isoscalar and isovector channels. The results concern calculations with T = pn pairs 
only and with all possible pairs taken into account. 
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